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Abstract 

A novel design of a state estimator is presented using 
second-order dynamic equations of mechanical systems. The 
eigenvalues and eigenvectors of the state estimator are as- 
signed by solving the second-order eigenvalue probletn of 
the structural system. Three design methods for the state 
estimator are given in this paper. The first design method 
uses collocated sensors to measure the desired signals and 
their derivatives. The second design method uses prefilters 
to shift signal phases to obtain estimates of the signal deriva- 
tives. These two methods are used to build a second-order 
state estimator model. The third design method is the con- 
ventional one which converts a typical second-order dynamic 
model to a first-order model, and then builds a state estimator 
based on the first-order model. It is shown that all the three 
deigns for state estimation are similar. A numerical example 
representing a large space structure is given for illustration 
of the design methods presented in this paper. 

Introduction 

Structural dynamic systems are generally described by 
second-order differential equations with symmetric and 
sparse structural matrices. Structural engineers and analysts 
perform dynamic analyses by taking full advantage of the 
symmetry and sparsity of the structural matrices to minimize 
the computational burden and keep physical insight intact 
For example, it is obviously easier to solve the eigenvalues 
of a symmetric and sparse matrix than a general matrix. On 
the other hand, control theory including estimation theory 
are established using first-order dynamic equations. Exist- 
ing control software tools today are thus written in first- 
order forms. In applications to structural dynamic systems, 
a composite state vector is used to transform the second- 
order dynamic equations to a larger dimensional first-order 
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form. Transformation to the first-order form not only in- 
creases the dimension by a factor of two but also destroys 
the symmetry of the structural matrices. As a result, signif- 
icant model reduction is generally required before any con- 
troller or. state estimator design can be accomplished, be- 
cause of the numerical difficulty associated with the solution 
of high dimensional equations such as Riccati equations. A 
number of researchers 1-8 have investigated model reduc- 
tion to circumvent the dimension problem. An alternative 
approach to model reduction is to preserve the second-order 
dynamic equations in designing the controller or state esti- 
mator. Recently, several researchers 6-10 have addressed the 
computational advantages of designing controllers and state 
estimators directly using the second-order structural models. 

The state estimator plays a major role in controller de- 
signs using state feedback under the constraint that the num- 
ber of sensors is less than the number of states. Second-order 
state estimator models have just recently received attention in 
the literature. An optimal state estimator known as Kalman 
Filter has been used in Ref. 6 for discretized second-order 
structural models. Robust computational procedures for solv- 
ing Kalman Filter estimation error covariance matrices have 
been developed for second-order models in Ref. 7. A dissi- 
pative state estimator in second-order form was introduced in 
Ref, 8. This state estimator was analogous to a dissipative 
controller® with collocated sensors and actuators whereby 
positive definite feedback gains were designed to insure sta- 
bility. The computational advantages of second-order state 
estimator models are discussed in Refs. 10 and 11. 

The. objective of this paper is to develop a robust state 
estimator for use with robust controllers. This research was 
stimulated by the work in Refs. 10 and 1 1 where sub-optimal 
second-order observers were developed using optimal control 
theory. The approach of this paper is to extend the technique 
presented in Ref. 12 for robust eigensystem assignment of 
second-order controllers. As in Ref. 12, the technique takes 
advantage of a second-order form of the system equations 
(instead of transforming to a first-order form) which results 
in coasiderable computational efficiency. The technique can 
handle any forms of feedbacks, i.e., displacement, velocity, 
and acceleration. It is known that the controller and state 
estimator in the first-order form are dual in a mathematical 
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sense, which implies that the design freedoms are identical. 
The question arises whether the same statement is true for the 
second-order form. It will be shown that the design freedoms 
for the second-order state estimators are only half of those 
for second-order controllers (see Ref. 12). Therefore, this 
paper presents methods of obtaining the additional freedoms 
which are needed to complete the design process. 

In this paper, three methods for the design of state estima- 
tors having second-order models are presented using eigen- 
system assignment techniques. The first design method uses 
collocated sensors to measure signals and their correspond- 
ing derivatives to gain full freedom to build a second-order 
state estimator model. The technique used in this design 
assigns the eigenvalues such that the resulting closed-loop 
system is robust with respect to system parameter uncer- 
tainties. This is accomplished by requiring the closed-loop 
eigenvectors to be as close as possible to the column space of 
a well-conditioned matrix. The second design method uses 
a second-order prefilter design in place of collocated mea- 
surements of the signals and their derivatives as required in 
the first design method. Here, the prefilter is designed to 
shift the phase of the signals, thus replicating the effects of 
the signal derivatives in the first design. The third design 
method is the traditional state estimator design in which the 
estimator is constructed based on a first-order model of the 
system. However, the gain matrix is computed through a 
second-order model. It is demonstrated that a second-order 
state estimator together with a prefilter design has the same 
design freedoms as the conventional first-order state esti- 
mator. A numerical example is given to demonstrate the 
proposed method. 

State Estimators with Measurement Signals 
and their Derivatives 

In vibration control of flexible structures, two set of 
second-order linear, constant coefficient, ordinary differential 
equations are frequently used. These equations, in matrix 
form, are 

Mb + Di + Kx= Bu (1) 

y=H v x + H d x. (2) 

Equation (1) is the system dynamic equation having x as the 
state vector of dimension n, and Af, Z?, and A" as the mass, 
damping and stiffness matrices, respectively, which generally 
are symmetric and sparse. The n x p influence matrix B 
describes the control force distribution for the p x 1 control 
force vector u. Equation (2) is the measurement equation 
having y as the measurement vector of length m, H v the 
m X n velocity influence matrix and B d the m x n deflection 
influence matrix. 

If the measurement vector y in Eq. (2) is used directly for 
a feedback control design, an output feedback controller is 
obtained. The output feedback control is generally attractive 
because it is simple and easy for real time implementation. 


However, a stable and robust output feedback controller 
may require either too many measurements which are not 
practical, or some measurement devices which are not yet 
available and need to be developed. On the other hand, 
the state feedback control law assumes that all states are 
measurable. In many practical control designs for flexible 
structures, it is physically or economically impractical to 
install the sensors that would be necessary to measure all the 
states. For such cases, a state estimator is needed to estimate 
the states from the measurement outputs, and provide enough 
freedom for a stable and robust feedback controller design. 

The basic approach of estimating the states is to simulate 
the state and output measurement equations of the system 
on a computer with an assumed initial state vector. In other 
words, Eqs. (1) and (2) are simulated on a computer with 
the same input u as applied to the actual physical system. 
For noise-free and uncertainty-free cases, the slates of the 
simulated system, i.e. the estimated states, will then be 
identical to the states of the actual systems if initial states 
are the same. However, the actual system may be subjected 
to unmeasurable disturbances which can not be used in the 
simulation but affect the output measurements. In order to 
make sure that the estimated state does not deviate too much 
from the actual state values, the difference between the actual 
output and the estimated output should be used as one of the 
driving inputs in the estimation equation. 

Let the state estimation and the output equations be 

Mb + Db + Kb = Bu + L d (y - y) + L v {y - y) (3) 

£ = H v b + H d b (4) 

where x is the estimated state vector of length n, y the 
estimated output of length m, and L d and L v the n x m state 
estimator gain matrices. Here, in contrast to the conventional 
approach, an additional term, Zt,(y-#), is added to Eq. (3) to 
penalize the difference between the actual output derivative 
and estimated output derivative. Why the additional term 
is added in Eq. (3) will be explained in detail later. Note 
that it is not wise to differentiate measurement signals to 
generate this additional term for real time implementation. 
Additional sensors are recommended for use in measuring 
the derivatives of measurement signals. 

To determine the matrices L v and L d , define the esti- 
mation error as 

e-x-b, (5) 

Subtracting Eq. (3) from Eq. (1) and employing the relation- 
ship given in Eqs. (2) and (4), the error equation becomes 

[M + L v H v ]e+ [J9+ L v H d + L d H v ]e + [K + L d H d ]e = 0. 

( 6 ) 

A question arises as to if an appropriate choice of the gain 
matrices Z* and L d will move the eigenvalues of Eq. (6) 
to the left-hand plane so that the steady-state value of e(t) 
for any initial condition is zero, i.e. lim e(<) = 0. The 
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following paragraphs present a novel way of synthesizing 
the matrices L„ and L d . 

Assume that the system, Eqs. (1) and (2), is observable. 
The left eigenvectors <& and eigenvalues A* for the system 
in Eq. (6) are related by the equation 


Jl {(M + L v Hv]*l + \D + L v H d + L d H v ] X k + [K + L d H d \ } = 0; 

k= l,...,n. (7) 

The subscript k refers to the mode number. Rearranging 
Eq. (7) in a compact matrix form yields 


[^k ^k (^k-^v + L d ) ] 


(X\M + X k D + K) 
(A k H v + H d ) 


= 0 ; 


k = l,...,n. 

The transpose of this equation gives 


( 8 ) 


tk 

{WT + Ll)*k 


— r k 4> k — 0; k — 1 , ...,n. (9) 

If the closed-loop eigenvalues X k (k - 1, ... ,n) and their 
complex conjugates are assigned, Eq. (9) can be used to 
determine the gain matrices L% , and L d . Because the vector 
<f> k is in the null space of the matrix T k , it is necessary to 
compute the null spaces of the matrices r* (k = 1, . . . , n) 
corresponding to the desired eigenvalues A* (k = 1 , . . . ,n). 
Tb obtain the nontrivial solution space of the homogeneous 
equation (9), the singular value decomposition (SVD) is 
applied to the matrix r* yielding 


r* = u k z k v: = u k 


Ok 

0 



( 10 ) 


Because Ak in T k is a complex value, all the quantities are 
complex except the diagonal matrix a k which contains the 
nonzero and positive singular values. Here the superscript * 
mm; transpose and complex conjugate. It follows that 
the matrix V ok represents a set of orthogonal basis vectors 
spanning the null space of the matrix r* so that 


Tk^k = r kKkCk = 0 01) 

where c* is an arbitrary column vector with an appropriate 
length. Note that if I* is well-conditioned (i.e. not close to a 
matrix of lesser rank which is easily found from the singular 
values; hence the advantage of using SVD), the above basis 
for null space V ok can be computed more efficiently by taking 
the QR decomposition of T*. If the matrix [X\M+X k D+K] 
is invertible, the vector $£ ^ k }*> where <f> k is an arbitrary 


vector ol length pand^ = —[X^M+X k D+K] *[A k Hj + 
Hj]<l> k , is in the null space of the matrix r*. 

To obtain an expression for gain matrices L v and L d , 
choose a particular set of vectors, 4> k = Vot c k (k = 
1,. . . ,n) satisfying Eq. (9), corresponding to some choice 
c k , and partition the vector <f> k into two components such 

r »4 s r, [^‘] s r, j^.°* J c, = o. 02 ) 

Comparison of Eqs. (9) and (12) yields 

A + = 4> (13) 


where the n x n matrix $ = (^ 1 ,^ 2 , . • • the m x n 
matrix 4 = • • • )$n] and the n X n matrix A = 

diag[Xi ,A 2 ,...,A„]. To solve for the gain matrices L v and 
L d - (Both matrices are real), decompose Eq. (13) into real 
and imaginary parts to yield 


Ar - *<Aj] + Ll<Sir = (14a) 

Ll (5 r A,- + Ar] + Lfii = (14b) 

or in matrix form 


L T H = [Ll L\ ) 


($, Ai + $,• Ar) (4>r Ar_J A,) ' 

<I>,' 


= [4< 4> r ] = - 


(15) 


Here, the subscripts r and i respectively refer to the real 
and imaginary parts of the associated quantities. The gain 
matrices L v and L d can be obtained from Eq. (15) 


L t = [L t v = 

. . , r(^A,- + $<Ar) ($ r Ar -*<A<)1“ 1 

= *r)\ *. * r 


L (16) 
A matrix inversion is required in the computation of the 
gain matrices Lv and However, if the number of eigen- 
values to be assigned is less than the number assignable, n, 
Eq. (11) becomes underdetermined which leads naturally to 
a minimum gain solution. To assure that the above matrix 
is well-conditioned for inversion, the condition number of 
the matrix ¥ should be the smallest possible. Interestingly, 
the above numerical requirement for the well-conditioning 
of the matrix inversion problem corresponds exactly to the 
eigenvalue conditioning problem since ® consists of eigen- 
values and eigenvectors. For sufficiently small damping 
D and small real part of the closed eigenvalue A rt in 
Eq. (15) approaches zero, because all the null spaces cor- 
responding to r* (k = 1, ... , n) are nearly in the real do- 
main. The_matrixjl' in Eq. (15) can be approximated by 
¥ = diag[$ r Ai $ r ]* In structural dynamics terminology, 
it indicates that for small damping and small real part of the 
eigenvalues assigned for the closed-loop model, the real part 
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of the closed-loop eigenvector matrix $ r dominates the con- 
ditioning of the matrix 'I'. The closed-loop eigenvectors are 
chosen and discussed in the following section. 

Eigenvector Assignment for Robustness 

Define an n x n well-conditioned matrix Ho, with vectors 
hoi,ho 2 , - ,tk>n as its columns. Then, the closed-loop 
eigenvectors are chosen to be as close as possible to the 
range space of the columns of matrix Hq to achieve a robust 
closed-loop design. If the open-loop conditioning is good to 
begin with, then the columns of matrix Ho may be chosen 
to correspond to the open-loop eigenvectors, i.e., 


hot = Ak ; * = 1,2, (17) 


In general, the choice of the open-loop eigenvector ip ok is 
arbitrary as long as the resulting closed^ loop , eigen vectors 
are linearly independent. If the control system is used 
only to provide active damping, and the closed-loop damped 
frequencies are quite close to their open-loop values, then 
it is best to use the open-loop eigenvector corresponding to 
an eigenvalue with the same (or similar) frequency for the 
vector V’ot in Eq- (17). 

Alternatively, the matrix Ho may be chosen to be an 
arbitrary unitary matrix (with perfect conditioning), or the 
closest unitary matrix to the open-loop eigenvectors. In 
the latter case, the matrix Hq is then the solution of the 
constrained least square problem minimizing 

|*o -#o| 

subject to HqHq = I 

which leads to 

H 0 = UW* (18) 


where U and W are, respectively, the left and right singular 
vectors of the open-loop eigenvector matrix ¥<» and I is 
an n x n identity matrix. It is noted that by choosing the 
closed-loop eigenvectors to be as close as possible to the 
column space of matrix Hq the closed-loop conditioning 
of the second-order system represented by Eq. (9) will be 
enhanced. However, an additional requirement is needed to 
ensure the well-conditioning of the actual system of Eq. (6), 
and that is to require that the estimator gains be as small as 
possible. 

Having defined Hq, the closed-loop eigenvectors 
^> k = V ok c k , k = 1 , 2 ,..., n, and the corresponding 
coefficient vectors c k are computed through the following 
sequential three steps. 


Step (1): Obtain the vector in the attainable closed-loop 
eigenvector space, (sec Eq. (12)), corresponding to the 
first closed-loop eigenvalue, which is as close as possible 

to the range space of the columns of matrix Qo = jjj® , 

where No is a p x n null matrix. The vector can be obtained 


from the algorithm described by Golub and Van Loan 13 for 
the computation of principal angles and vectors of a subspace 
pair. Expand both matrices T 0 i and Qo in terms of their QR 
decomposition, i.e., 

'Poi = 0v o i^v cl > Qo = Qq 0 Rqo' (1®) 

Here Qy and Qq 0 are orthonormal matrices of dimensions 
(n+p) xpand (n+p) xn, respectively, and Hy ol and Rq 0 
are p x p and n x n upper triangular matrices, respectively. 
Project the vectors Qy ol unto Qq 0 to obtain 

H = Q' Vcl Q Qo . (20) 

The singular values of matrix H are cosines of the principal 
angles of the subspace pair {^(Q^^^Qqo)}. 15 Where 
!ft( ) denotes the range space of ( ). Taking the singular 
value decomposition of H, gives 

H = YSZ* (21) 


in which Y and Z are unitary matrices, and the matrix 5 
contains the singular values of H, i.e.. 



Note that the singular values o\, . . . , o v are all positive and 
less than or equal to 1 because the matrix H is formed by 
orthonormal matrices as shown in Eq. (20). The angles 
0j,. ..,6 P are -file principal angles of the subspace pair. 
Having in mind that 1 > oi > 02 >...> cr p , the vector 
in the orthonormal columns of Qy 0l which is closest to the 
column space of Qq 0 is then 

T>x = Qv 0 ,y = VoxR=r\v (23) 

in which y represents the first column of matrix Y . Choosing 
the first closed-loop eigenvector to be the coefficients ci 
of Eq. (12) become 

Cl = ifc 1 y. (24) 

Y ©1 

Note that Jfcr defined in Eq. (19) is a nonsingular matrix. 

Step (2): Reduce the column space of matrix Qo by the 
vector <^i, in order to ensure the linear independence of the 
closed-loop eigenvectors. If ^ is not in the column space 
of Qo, then reduce the column space by the closest vector 
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in the space to the vector The new subspace is spanned 
by the columns of matrix Q\ defined as 

Qi = [Qo-Qq 0 ^Qq,Qo) <25) 

where z is the first column of matrix Z. Note that the 
columns of matrix Qi are orthogonal to the vector Qq 0 z, i.e., 
z*Qq o Qi = 0. Equation (25) may not be computationally 
efficient It is shown here just for simplicity and clarity, and 
other computational procedures to compute Q\ may be used 
instead. 

Step (3): Repeat steps (1) and (2) for the remaining n- 1 
eigenvalues. For the ? h eigenvalue the column space of 
matrix Qi-i is reduced by the vector Qq { _ x z resulting in a 
new matrix Qi, 

Qi = [Qi - 1 - QQ i - l **'Qq i - l Qi-& ( 26 ) 

In summary, following the developments outlined in steps 
(l)-(3), the closed-loop eigenvectors are chosen to be as close 
as possible to the column space of a well-conditioned matrix 
Hq and the estimator gain matrix elements are designed to be 
as small as possible, thereby resulting in a robust closed-loop 
design. 

Note that the above formulations are nearly identical to 
the eigensystem assignment with full state feedback 12 except 
for slight differences in the matrix I'jt defined in Eq. (9), 
i.e. the matrix B in the case of full state feedback is 
replaced by matrix [X k Hj + Hj] for the state estimator. 
Computationally, both state feedback and state estimator 
designs are identical. 

Let us come back to discuss the term Lo(y — ij) which 
was yfldftd in Eq. (3) to gain more freedoms and make the 
estimator problem dual to the state feedback problem. Exam- 
ination of Eq. (13) reveals that when 1% = 0, Zj = 1 

which is, in general, a complex matrix. It thus contradicts 
the requirement that the gain matrix Zj must be real. Note 

that is function of the closed-loop eigenvalues of the 
state estimator. Consequently, it is immediately concluded 
that, with the absence of Zj', the solution of the gain matrix 
Z in real domain, in general, does not exist However, there 
may exist certain eigenvalues for the state estimator such that 
4<F -1 is real. For example, given a set of desired estimator 
eigenvalues, the gain matrices, Lv and L d , are computed 
from Eq. (16). Using the approach shown in Ref. 10, L v 
is omitted from Eq. (6), and the eigenvalues of the system 
in Eq. (6) may be still in the left-half plane. However, the 
eigenvalues thus obtained are, in general, different from the 
desired ones which are originally used to compute the ma- 
trices Z„ and L d . However, if the eigenvalues satisfy the 
performance requirements, it can be used for real time im- 
plementation. To this end, it is concluded that, using the 
second-order dynamic model for design of state estimators. 


requires an additional expense in the sense that collocated 
sensors may be needed for measuring the signals and their 
derivatives. 

State Estimators with Measurement Signals 
and Prefilters 

In some cases, collocated sensors for measuring signals 
and their derivatives are not available. An alternate way 
is needed to gain the freedom for designing the state esti- 
mators using second-order dynamic models. The purpose 
of adding sensors, as discussed previously, is to measure 
signal derivatives to obtain signals with different phases so 
as to gain enough freedom for the state estimator design. 
There are many other ways to shift signal phases without 
additional sensors. An alternative is the use of prefilter to 
approximately estimate the signal derivatives. 

Let- the state estimation and the output equations be 

Mi Dx + K± = Bu + L d (y -y) + L v (z - z) (27) 

z + Pz + Qz = y (i.e. H v x + H d x) (28) 
and 

z + Pz+Qz = y (i.e. H v i + H d i) (29) 

where the vector z of length m contains the filtered signals, 
z the estimated output vector of length m, and Z„ the n x m 
gain matrix associated with the error between z and z. Here, 
the additional term, L v (z-z ) provides the freedom necessary 
to design the state estimators for second-order models. The 
phase differences between the measured signals y and the 
shifted signals z are determined by the mxm square matrices 
P and Q. For simplicity in real time implementation, the 
matrices P and Q may be chosen to be diagonal such that 
signals are not coupled in the prefiltering processes. Indeed, 
in this case, a scalar second-order equation is obtained for 
each signal which can be easily implemented by an analog 
computer. Subtracting Eq. (27) from Eq. (1) with the aid of 
Eqs. (28) and (29) yields the error equation 

Me + [D + L d H v )e + [K + L d H d ]e + U e, = 0 (30) 


and 

'e t + Pe K + Qe t -H v e-H d e = 0 (31) 

where e M = z - & is the error between the filtered-signals 
and the estimated filtered-signals. Equations (30) and (31) 
can be rewritten in the following matrix form 


M 

0 


01 /e\ D+L d H v 

JjVeJ + [ ~H V 



+ \ K -*‘ q] (e!)=°- < 32 > 

The left eigenvectors 4%] (k = 1, . . . , n + m) and the 
corresponding eigenvalues A* for the system of Eq. (32) are 
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related through the following equation 

r m o 


of [ o /] x i 


D+L d H v 0 
-H v P 


K + L d H d Lv 1 "I _ 
-H d Q\j-° 


{[? !]*♦ 


Ll 


’D+HflT -Hf 
0 P T 


= 0. 


-HJ] 

11 

r&i 

Q T 

II 

h' 


(33) 


where subscript k refers to the eigenvalue number. The 
transpose of this equation is 


(34) 


This equation can be decomposed into two parts which are 


design parameters for the prefilter equations, they may be 
chosen to be diagonal for simplicity in real time implemen- 
tation as well as computation of L v and L d in Eq. (37). 


First-Order State Estimator Models 

There are two reasons why first-order state estimators for 
flexible structures are presented in this section. First, a com- 
putational procedure which is different from conventional 
ones is developed hoe. Second, a comparison between the 
first-order and the second-order state estimators is given for 
better understanding of the characteristics and merits of both 
approaches. 

Equation (1) can be rewritten in a first-order form, 


7 0 ' 
0 M 



(38) 


{xlM+XtiD+HtLft + lK+HjLtyh 
-[A t J?J + J?J]& = 0 

and 

4’3k + [Ai + A*P r + Q r ]& = 0. (35b) 

Solving for from Eq. (35b) and substituting it into 
Eq. (35a) yields 


[(XlM+X k D+K T ) (A k H? + Hj)) 
\ T>k 

•|{iJ + [Ai + A k P T + Q T ]- l Ll)t k 


= ^kfa = 0 . 

(36) 


This equation is nearly identical to Eq. (9) in the sense 
that iys in both equations are identical, and therefore 
have the same null space. This simply means that the 
computational procedure developed previously can also be 
used in computing the gain matrices L v and L d . However, 
the way to compute L v and L d in Eq. (36) is somewhat 
different from that of Eq. (9). 


If the closed-loop eigenvalues A* (k = 1, . . . ,n + m) 
and their complex conjugates are assigned, the gain matrices 
can be determined as follows. As shown in Eq. (12), choose 
a particular set of vectors fa (k= 1, . . . , n + m) satisfying 
Eq. (36) and partition vector fa into two components, then 
Eq. (36) implies 


{L d 4- [Aj[ -f AfcP T + Q t ] 'Ll3k — $k 
k = 1 ,...,n + m. 


(37) 


This equation can be decomposed into real part and imagi- 
nary parts to solve for the gain matrices L v and L d (Both 
matrices are real), containing 2m x n unknown elements. 
However, the computational procedure is not as straightfor- 
ward as the previous case, Eq. (13). Because P and Q are 


Correspondingly, the state estimation equation becomes 


°1 
0 M 



[-* -9M 
[£][* *](!:*) 


u 

■ (39) 


Note that the gain matrices L v and L d used here are 
somewhat different from those shown in the above section 
(Eq. (3)), even though they look similar. To determine the 
matrices L v and L d , the state estimation error is defined by 


e T = [(x - &) T { x — x) r ] r . (40) 


Observe that the estimation error defined in Eq. (5) for a 
second-order model is different from the one defined above, 
Eq. (40). There is no estimation error for velocity terms 
involved in Eq. (5). Instead, the estimated velocity errors 
are incorporated in the state estimation equation, Eq. (3). 
Subtracting Eq. (39) from Eq. (38) yields the error equation 


■M 0 
0 I 


K D e+ l Sd ^*'] e = 0 -^ 41 ) 


The left eigenvectors {4% u 4$ d ] {k = l,...,2n) and 
the corresponding eigenvalues A* for the system given by 
Eq. (41) are written as 



where subscript k refers to the eigenvalue number. The 
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transpose of this equation is 



This equation can be decomposed into two parts which are 


*k*kv + K*kd + *d \$*kv + L&kd] = 0 («*) 

and 

" tkv + Dfkd + Hv [f'vVfcv + tkd] = (***>) 

Premultiplying Eq. (44b) by A* and adding the resulting 
equation to Eq. (44a) yields 


[a£m + A k D + if] 4kd + [*k H I + H d] { L v*kv + Ld+kd] = 0 

or in a compact matrix form 


[(XlM+X k D + K) (A k H* + fif)] 

f 4>kd 

'\(Ll<t> kv + LT<f> kd ) 


z=T k <j> k = 0. 


(46) 


This equation is nearly identical to Eq. (9) in the sense that 
iys in both equations are identical, and therefore have the 
same null space. This is a significant result which simply 
inHimtp-c that the computational procedure developed above 
for the second-order model can also be used in computing 
the gain matrices L v and L d for the first-order model. The 
matrices Lv and L d can be determined when all the null 
spaces of the matrix T k corresponding to the eigenvalues A* 
are computed. However, the way to compute L v and L d in 
Eq. (46) is somewhat different from those in Eq. (9). 

If the closed-loop eigenvalues A* (k = 1 ,. . . ,n) are as- 
signed including their complex conjugates, the gain matrices 
Lv and Xj can be determined as follows. Following Eq. (12), 
choose a particular set of vectors 4> k (k = 1, . . . ,n) satisfy- 
ing Eq. (46) and partition the vector 4>k into two components, 

<tf[ = Equation (46) implies 

L[ 4>kv + kd = $k i k = 1 , . . . ) n (47) 

where <f> k „ can be solved from Eq. (44) 

<f> kv — — (K4> k d + -Sj^*)/Afc '> k= l,...,n. (48) 

Both Eqs. (47) and (48) arenxn equations compared to a 
2 n x 2n equations solved for a typical first-order model such 
as Eq. (43). In other words, there are 2nx n less equations to 
solve for the matrices Lv and L d for each eigenvalue using 


Eqs.(47) and (48). Therefore, a total of 2n x n x n equations 
are computationally saved for n assigned eigenvalues. _Now 
let $ = [<Aidi- • • >^wl* $ = • • • >^»») and ^ = 

[& , . . • , $„] • Equation (47) gives 


Xj* r + X^ r = 6 r (49a) 


I**, + IR- = (49b) 


or in matrix form 


m [!; I:] = i*< 


where the subscript r and i respectively refer to the real 
and imaginary parts of the associated quantities. The gain 
matrices X„ and L d can then be solved by using 


[Ll XJ]=[^ *r] 


j *< 



(51) 


Again, a matrix inversion is required in the computation of 
the gain matrices L v and L d . 


Numerical Example 

The second-order state estimator is used to control the 
vibrational modem of the flexible truss structure shown in 
Fig. 1. The structure genetically represents a test article for 
the NASA’s Controls-Structures-Interacdon program. It is 
composed of a L-shaped bus, a reflector, a laser feed, and 
two suspension cables used to simulate on-orbit conditions. 
The original finite element model is composed of 350 grid 
points with six degrees of freedom per grid point, resulting 
in a 2100 degree-of- freedom. However, for preliminary 
control designs a reduced order model comprised of the 
first nine modes of the structure, covering a bandwidth of 
0-bBz, is used. These modes include four pendulum modes 
(modes due to suspension effects) and five flexural as well 
as torsional modes. The actuation of the control forces is 
provided through six proportional gas thrusters located at 
various points on the body as illustrated in Figure 1. Six 
inertial measurement units and six accelerometers provide 
twelve measurements of linear velocity and accelerations 
at six locations along the body (almost collocated with the 
actuators) as indicated in Fig. 2. 

The open-loop eigenvalues and the desired closed-loop 
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eigenvalues are summarized in Table 1. 


Thble 1. 

Open-Loop Closed-Loop 

Eigenvalues Eigenvalues 


0.0000 + 0.7746t 
0.0000 + 0.8019* 
0.0000 + 0.8043t 
0.0000 + 3.9237* 
0.0000 + 13.4800t 
0.0000 + 15.2194* 
0.0000 -(- 18.9833t 
0.0000 + 19.4861* 
0.0000 + 23.4345t 


-0.0501 + 1 .0000* 
— 0.0601 + 1.2000i 
-0.0751 + 1.50002 
-0.1964 + 3.9237* 
-0.6748 + 13.4800* 
-0.7619 + 15.2194* 
-0.9504 + 18.9833* 
-0.9755 + 19.4861* 
-1.1732 + 23.4345* 


The closed-loop eigenvalues are chosen such that the first 
four pendulum modes are provided with 30% damping and 
the remaining modes with 5% damping. The closed-loop 
damped frequencies (imaginary part of the eigenvalues) are 
larger than the open-loop values for the first three modes, but 
the same for the remaining modes. The closed-loop eigen- 
values are assigned via a constant full state feedback, using 
the robust second-order assignment technique described in 
Ref. 12 in conjunction with the first design method described 
in the section of State Estimators with Measurement Signals 
and their Derivatives. Here, however, the estimated state is 
used instead of the actual state in the feedback loop. 

Using the closest unitary matrix to the open-loop eigen- 
vector matrix as the choice for matrix Hq, the eigensystem 
assignment technique results in a well-conditioned closed- 
loop system with a modal matrix condition number, c(^ ), of 
26.76 which compares quite well with open-loop condition 
number of 23.43. The forbenius norm of the gain matrix 
is quite small at 6.89. These results indicate that proposed 
eigensystem assignment technique is quite effective and can 
lead to viable closed-loop designs. 

The estimator gains are also obtained using the eigen- 
system assignment technique described in the first design 
method. The closed-loop eigenvalues of the estimator sys- 
tem are chosen equal to the closed-loop eigenvalues of the 
actual system given in Table 1 except that the real part of all 
the eigenvalues are chosen at —0.5 to achieve an acceptable 
performance for the overall system. The resulting state esti- 
mator has a good conditioning of 67.06 and a low norm of the 
gain matrix at 1.10 which further illustrates the effectiveness 
of the proposed design procedure to obtain well-conditioned 
closed-loop systems with small control effort. 

To verify the feasibility of the designed state estimator, 
a numerical simulation was carried out wherein the dynamic 
behavior of the closed-loop system for an initial disturbance 
is investigated. The time history for the first pendulum mode 
and its associated estimation error are presented in Figs. 2(a) 
and 2(b), respectively, for an intial velocity of 1.0 in all 
coordinates. Similarly, the time histories for the first flexible 


mode, and its corresponding estimation error are given in 
Figs. 2(c) and 2(d). The time histories of the control forces 
are respectively illustrated in Figs. 3(a)-3(f). The results 
indicate that the initial disturbance is practically damped out 
in 10 — 12 seconds. 

Concluding Remarks 

Three design methods for a state estimator were presented 
in this paper. The first two methods were used to build a 
second-order state estimator model. The third design method 
was the traditional state estimator design using a first-order 
model, but the gain matrix was computed through a second- 
order model. Careful examination of the the third method 
reveals that the first-order state estimation equation does in- 
clude a filter equation. Indeed, when a second-order model 
is converted into a first-order model, an additional first-order 
equation is generated which is then implicitely used to build 
a filter equation. Consequently, it can be concluded that 
as long as a prefilter design is added in the second-order 
dynamic model, frill freedom to design a state estimator is 
obtained. From the computational point of view, the second- 
order models are more attractive for use in designing the state 
estimators, because the dimension of the mathematical mod- 
els remain unchanged, rather than an increase by a factor of 
two for the first-order models. Furthermore, the fundamental 
structure of the mathematical models such as the symmetry 
and sparsity of the mass, damping and stiffness matrices is 
maintained. The disadvantages of the state estimation using 
second-order models include the requirements of additional 
sensors or prefilters for real time implementation. 
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